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Abstract 

An atmospheric-correction method which uses cloud-shaded pixels together with 
pixels in a neighboring region of similar optical properties is described. This cloud-shadow 
method uses the difference between the total radiance values observed at the sensor for 
these two regions, thus removing the nearly identical atmospheric radiance contributions to 
the two signals (e.g. path radiance and Fresnel-reflected skylight). What remains is largely 
due to solar photons backscattered from beneath the sea to dominate the residual signal. 
Normalization by the direct solar irradiance reaching the sea surface and correction for 
some second-order effects provides the remote-sensing reflectance of the ocean at the 
location of the neighbor region, providing a known "ground target" spectrum for use in 
testing the calibration of the sensor. 

A similar approach may be useful for land targets if horizontal homogeneity of scene 
reflectance exists about the shadow. Monte Carlo calculations have been used to correct 
for adjacency effects and to estimate the differences in the skylight reaching the shadowed 
and neighbor pixels. 

Key words: Satellite calibration, cloud shadow. 
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Introduction 


Accurate calibration of aircraft and space-borne sensors that view dark targets such as 
the ocean is critical. As much as 90% of the signal at the sensor is due typically to the 
atmosphere; as little as 1 0% of the signal may be due to target reflectance 1 - 2 . Therefore, a 
5% sensor-calibration error may result in an error of as much as 50% in the calculated 
reflectance. Prelaunch calibration accuracies of 2-5% are representative of requirements 
for space sensors, and in-orbit sensor performance typically differs from prelaunch 
performance. Since sensor performance degrades over the life of the mission, occasional 
recalibration using aircraft- or ground-based methods 1 - 4 is required. 

The method of choice for sensor calibration when viewing the ocean from aircraft or 
space is to locate a region with relatively stable and homogeneous optical properties (e.g. 
the Sargasso Sea), measure its water-leaving radiance spectrum, and determine the optical 
properties of the atmosphere coincidentally with a sensor overpass 4 - 2 - 5 - For coastal 
applications with variable bottom reflection and horizontal gradients in water constituents, 
locating a region where surface calibration can be accurately performed is often quite 
difficult. 

In this work we describe an atmospheric-correction method that uses cloud-shadowed 
pixels in combination with unshadowed pixels in a neighboring region of similar optical 
properties. The cloud-shadow method uses the difference between water-leaving radiance 
values L W (X) for these two regions. This allows nearly identical contributions to the two 

signals (e.g. path radiance and Fresnel-reflected skylight) to be removed, leaving mostly 
solar photons backscattered from beneath the sea to dominate the residual signal. 
Normalization by the direct solar irradiance reaching the sea surface provides, to first order, 
the remote-sensing reflectance of the ocean at the location of the neighbor region. Special 



attention must be paid to evaluating and correcting for adjacency effects and the difference 
in skylight reaching the shadowed and neighbor pixels. 

The Airborne Visible-InffaRed Imaging Spectrometer (AVTRJS), flying at an altitude 
of 20 km, provides data simulating that expected from hyperspectral space sensors of the 
future. 6 Data from AVIRIS were used to develop and test the methodology described in 
this paper. 

Theoretical Considerations 

For illustrative purposes, imagine a viewing scenario in which the sensor calibration is 
correct. The solar zenith angle, 0 O , is 45°, and the angle from the pixel to the sensor, 0, is 

about 0° (see Fig. 1 ). A small, compact cumulus cloud removes direct solar photons and 
shadows a region. The water-leaving radiance directed toward the sensor from this 
shadowed region is designated L ws . (Note that terms indicating wavelength-dependence 

have been left out for brevity in cases where doing so is unlikely to cause confusion.) This 
radiance results from skylight photons reflected by the surface or scattered from beneath 
the ocean surface. 

Adjacent to the shadowed region is a neighboring patch of water with inherent optical 
properties identical to those of the shadowed region. This region is illuminated by direct, 
solar photons as well as skylight . The water-leaving radiance from the neighbor region is 
designated L^. 

In addition to the water-leaving radiance, the sensor measurement includes the effects 
of path radiance caused by atmospheric scattering of photons into the field of view of the 
sensor. The primary processes responsible for path radiance are molecular, or Rayleigh, 
scattering, and particulate, or aerosol scattering. Path radiance may be attributed, therefore, 
to photons which have suffered only Rayleigh scattering, only aerosol scattering, or some 
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combination of both. These three types of path radiance are denoted by L r L a , and L ra , 
respectively. 

Accordingly, let the total radiance measured at the sensor when viewing a neighboring 
area in unshadowed water be given by 

^tn ~ L r + L a + L ra + tjL^ , ( 1 ) 

where t d represents the diffuse transmittance of the atmosphere for water-leaving radiance 7 . 

The radiance measured at the sensor when viewing a shadowed pixel may be 
expressed in the same form, but some differences in the path radiances and diffuse 
transmittance may be expected. We assume that the cloud is of sufficient thickness that the 
direct solar beam is completely occluded. Since part of the viewing path to the shadowed 
pixel is also shadowed, this portion of the viewing path must produce less path radiance 

(Fig- 1). 

The apparent path transmittance of the water-leaving radiance from the shadowed 
pixel may not be equivalent to the term used in Eq. (1). Use of the "diffuse transmittance" 
is justified when viewing a large, homogeneous area. In such a case, target radiance 
scattered out of the viewing path is balanced by the radiance scattered into it from adjacent 
areas of the scene. In the case of the shadowed pixel, the adjacent areas of the scene are 
generally brighter, and so the apparent transmittance of the viewing path to the shadow 
will be enhanced by photons emanating from the bright portion of the image and scattered 
into the field of view of the sensor. 

With these ideas in mind, the total radiance at the sensor when viewing a shadowed 
pixel may be written 

F ts — F r - AL r + L a - AL a + L ra - AL ra + ( t d + , (2) 
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where A terms express perturbations due to nonhomogeneity in the scene illumination. 


In general, water-leaving radiance is the result of backscattering of light which has 
penetrated the air-sea interface, and may be expressed as the sum of two parts: one part 
caused by backscattering of diffuse skylight, and the other by backscattering of the direct 
solar beam. For the neighbor and shadowed pixels, respectively. 


^wn sky^wn sol^wn 


^ws skv^v 


(3) 


since sol L ws = 0. 


Even though the cloud is small and occludes a similar portion of the sky for both the 
shadowed and neighbor pixels, the diffuse irradiances incident cm the two pixels are 
unequal (see appendix). The cloud occludes the brightest part of the sky from the 
shadowed pixel, i.e. the part containing the radiance due to near-forward scattering by 
aerosols. On the other hand, the neighbor pixel is illuminated by the relatively bright side 
of the cloud. Therefore, the diffuse irradiance at the neighbor pixel may be greater than 
that in the shadow. Under the assumption that the remote sensing reflectance for skylight 
is about the same at the shadow and neighbor pixels, we may write 


sky^wn sky^ws Asky^-'u 


(4) 


so that 

Ltn ‘ ^ts — + ^a + ^ra + (sol^wn + Asky^ws) " sky^ws ■ (3 ) 

The first three terms on the right of Eq. (5) depend cm the length of the shaded portion 
of the viewing path to the shadow pixels. The height of the intersection of the viewing 
path and the upper edge of the cylinder of atmosphere shaded by the cloud can be 
determined from scene geometry. The layer of atmosphere below this intersection would 
be the source of less than 1 5% of the Rayleigh scattering, so we will assume that the 



Rayleigh-aerosol correction term, AL ra , is negligible. Then, following Gordon et al. 8 , for 
the aerosol and Rayleigh corrections we may write 

AL x = l“x T 'x F "o p x ( 0 - 9 o)} 1 ’d /47C . x = a,r, (6) 


where 

t' x = optical thickness of shaded viewing path for process x; 

G> x = single-scattering albedo for process x; 

F 0 " = F 0 exp( -(x - t’)/cos( 0 o ) ) , direct beam incident at height of shadow/path 
intersection; 

F 0 = extra-terrestrial solar irradiance; 

P X (0,0 O ) = { P x (0.) + [p(0) + p(0 o )] P x (0+) } / cos(0); 

cos(0+_) = +-cos(0 o )cos(0) + sin(0 o )sin(0)cos(<jHj> o ); 

p(0) = Fresnel reflectance for incident angle 0; 

P x (0) = scattering phase function for process x; 

f d = exp{-( (x r - x' r )/2 + (x^ - t'^) )/cos(0) }, diffuse transmittance from shadow/path 
intersection to top of atmosphere. 

A term involving (x a - x'J and was neglected in the last expression because diffuse 

transmissivity due to aerosols from the top of this layer is approximately unity 2 . Also, 
since ozone is generally found above the cloud layer, the term x'^ is also negligible. 

The term involving At d in Eq. (2) represents the apparent increase in diffuse 
transmittance of water-leaving radiance when viewing a shadowed pixel. The work of 
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Tanre et al. 9 indicates that this term is proportional to the difference in water-leaving 
radiances from the shadowed and neighbor pixels, and is dependent on the geometry of the 
particular case. We will express At d as 

— ® ( ^wn “ ^ws ) ^ Lws — ^d ® ( sol^wn Asky^ws ) ^ sky^ws . ( 7 ) 

Unlike for the AL x (X) terms, the values for ct, and hence At*(X), can not be determined 
using Lowtran 7 10 . For this application, a was estimated using approximations 
developed in Reinersman et al. u (Fig. 2). These approximations are based on linear 

scaling of Elterman’s 12 atmosphere to match the prevailing atmospheric aerosol optical 
thickness. To get an initial estimate of t a (X), a scaling relationship between L a (780 nm) 

and the aerosol optical thickness, t a (780 nm), was developed using Lowtran 7 (Fig. 3). 
First, L t (780 nm) was calculated for several values of r a (780 nm). Then L,(780 nm) was 

calculated for an atmosphere characterized by Rayleigh scattering only. Subtraction of the 
latter value from the L t (780 nm) values calculated for the atmospheres containing aerosols 
yields L a (780 nm) because L r (X) + L ra (X) is approximately constant for x a (X) < .5 . Based 

on the above simulations, we estimate that a(550 nm) is about 0.06 for shadows of 
spherical clouds of the size modelled in this paper when viewed through very clear 
atmospheres ( t a (550 nm) = .0375 ). 

Using Eqs. (7) and (5) we obtain 

+ (L u -L 0 -AL,-AL r }/t d (l-o). (8) 

Now define 8 

e(Xi,Xj) = { o> a (X,) x' a (X,) P a (0,e o ,X,) } / { o) a (Xj) fjl,) PJOfioty } , (9) 

and 

S(K Xj> = AL.ft) / AL.jJ.,) - V { F"A) W } / { F" 0 (Aj) }. (10) 
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Using these results in Eq. (8) 


t d (W - <*(>•,)){ + * ky L w A) } = LJX.) - L ts (X.) - ALA) 

-S-(X ; , Xj){ LJXj) - L ts (Xj) - AL/Xj) - t d (X,)(l - a(X ,)K JLJfy + ^Aj)) > • 

( 11 ) 

If Xj is a wavelength such as 780 nm where the water-leaving radiance is essentially 
zero, then ^L^A,) + Ask y L w S (^j) = and Eq- (5) becomes 

AL a (780nm) = L tn (780 nm) - L ts (780 nm) - ALA80 nm) . (12) 

This allows Eq. (1 1) to be written 

tA)(l - oCMI + Asky L ws (X,) } = L ta (X i ) - LJX.) - ALA,) 

-S A> 780 nm){ L ta (780 nm) - L ts (780 nm) - AL,(780 nm) } (13) 

From scene geometry, the cloud height, h, and height of the shadowed viewing path 
can be calculated. Then, ALA) can determined for all X by subtracting Lowtran 7 

results for the air column down to height h from the results for the entire air column. The 
aerosol correction term, AL a (780 nm), is then determined by Eq. (12). Knowledge of wind 
speed allows estimation of e(X j; 780 nm) 13 , and thus S'(Xj,780 nm). Asky L w A) Tm Y t> e 
estimated by Monte Carlo simulation (see appendix). Thus, Eq. (13) allows ^jL^A) to be 
calculated. 

If we assume a Lambertian sky, the average cosine for irradiance at the surface due to 
skylight, sky E d (0 + ), is 0.707 14 . This is roughly equivalent to having all skylight photons 

striking the sea surface at 45°, and is the same as the cosine of the irradiance due to direct 
solar illumination, (0 + ) , in this example. So, since reflectance is independent of the 
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source intensity and color, and illumination geometry is equivalent for the solar and 
average sky photons, 

A = »iLwn / *,.E d (0 + ) = L w / E d (0 + ) = R.J, (14) 

where is the remote-sensing reflectance due to total downwelling irradiance, E d (0 + ). 
Thus, deriving in this manner provides a calibrated target reflectance value that can be 
used to derive the atmospheric aerosol characteristics in a manner similar to that used by 
Gordon et al. 8 for low-chlorophyll, offshore waters. 

Experiment 

The cloud-shadow reflectance method was tested using Airborne Visible-InfraRed 
Imaging Spectrometer(AVIRIS) data collected from an altitude of 65,000 ft over the Straits 
of Florida. AVIRIS produces images consisting of 5 1 2 rows of 614 cross-track samples. 

At each sample location radiance is measured simultaneously in 224 channels. Spatial 
resolution for nadir viewing is 20 m. 

We calibrated the Airborne Visible-InfraRed Imaging SpectrometerfAVTRIS) at a 
clear-water offshore location in the Florida Current using the vicarious calibration method 
of Carder et al. 5 In brief, the R^ curves for calibration scenes are measured; the values of 
E d (0 + ) are obtained from Lowtran 7; and L W (X) values are calculated by multiplying R^ by 
E d (0 + ). This technique maintains consistency between the illumination for the atmospheric 
correction program and the ocean measurements. The extra-terrestrial solar radiance F 0 

provided by Neckel and Labs 15 has been used as the solar source for all model calculations. 
Then Lowtran 7 was used to calculate the total radiance, L t , reaching the sensor from the 

ocean and the atmosphere. The resulting data were median-filtered using the brightness at 
780 nm to discriminate against some 10% of the pixels apparently containing white caps 
and/or sun glint. Crests of waves of about 1 00 m in wavelength were observed in the 
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imagery, and this enhanced brightness would have been misinterpreted by the program as 
being part of the atmospheric path radiance. Finally, the calibration of the sensor was 
adjusted so that modelled and measured L t values matched. The first test of the cloud- 
shadow method was verification that the correct R^fX) = L W (X) /E d (X) spectrum could be 
recovered from data acquired from an improperly calibrated sensor. L t (X) spectra from the 
bottom cloud-shadow pair of Fig. 4 (Cloud 1), calibrated as described above, were used to 
perform the verification. The test procedure was as follows: 

1 ) The "true" L t (>,) spectra from both the shadow and neighbor pixels were increased 
by 10% to simulate spectra, L' t (X), which would have been acquired had the sensor 
calibration been in error by 10%. 

2) OT ,L' WT1 (A.) was calculated from Eq. (5) as described earlier using L' t (X). Then the 
first estimate of the remote-sensing reflectance spectrum, R' re (A.), was obtained from 

= 05 ) 

3) Corrected L' t (X) spectra at the aircraft were simulated by adding atmospheric effects 
derived from Lowtran 7 for the appropriate conditions, i.e. , a nadir-viewing sensor, a 50.8° 
sun angle, a 60 km visibility through a marine aerosol with 80% relative humidity and 
1014 mb atmospheric pressure. 

The iteration process consisted of repeating steps 2 and 3 above. The results of this 
method are shown in Fig. 5. Both the "true" and "miscalibrated" reflectance spectra are 
illustrated along with the intermediate results from the cloud-shadow method. After four 
iterations of recalibration, the "corrected" cloud shadow reflectance spectrum approximates 
the "correct" spectrum as closely as if the correct calibration had been used initially. 

The method is effective for two reasons: 1 ) the error in the solar spectrum used to 
illuminate the atmosphere in the Lowtran 7 calculations is less than 2%; 2) the atmospheric 
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effects simulated by Lowtran 7 account for most of the radiance L, measured by the sensor, 

and for clear days this can be very accurately determined. If the sensor calibration were 
10% too large, the L w spectrum calculated conventionally from the difference between L t 
and (L a + L r ) would contain nearly all of the error. Thus, the conventionally calculated L w 

could be 40% or more higher than the true curve. The atmospheric path radiance removed 
inherently by the cloud-shadow method leaves an estimated L w spectrum that is high by 

only about 10%, however. Since this spectrum is closer to the correct one, recalibration 
based upon this new ground-target reflectance provides a better basis for a second iteration 
of the recalibration loop, providing a calibration factor that is in error by less than 5%. With 
iteration, convergence toward a calibration factor consistent with the solar spectrum and 
the model atmosphere used in Lowtran 7 is assured. 

Results and Discussion 

Two further demonstrations of the cloud-shadow atmospheric-correction method were 
performed using the scene shown in Fig. 4. The data were collected on a SE-NW transect 
from the Florida Current to Biscayne Bay. Elliot Key is shown in the upper left comer 
offshore. The shadow-neighbor pairs are near the seaward reef edge of Biscayne National 
Park. The cloud image associated with the bottom pair (Cloud 1) is recorded on the 
adjoining scene. Moments before this imagery was collected, AVERIS acquired data for 
calibration over the Florida Current from a scene some 25 km to the southeast. 

The results of applying the cloud-shadow method to Cloud 1 are compared to 
conventionally derived results in Fig. 6. Results shown for the conventional method were 
derived as described in Carder et al. 5 with one exception: the water-leaving radiance at 780 
nm could not be assumed to be zero, as the method usually requires, because of sunglint 
and/or foam effects due to the 10-12 m/s winds 16 . Instead, the water-leaving radiance at 
780 nm was assumed to equal that determined by the cloud-shadow method, since it 



inherently separates the atmospheric radiance from all water-leaving radiance including 
sunglint and foam effects. The reflectance derived using the cloud shadow method differs 
from the conventional result by less than 1 0% from 400 nm to about 460 nm. From 460 
nm up to 900 nm agreement is within 5% except for a minor perturbation at about 765 nm. 

The results of applying the method to the top shadow-neighbor pair (Cloud 2) are 
illustrated in Fig. 7. For this region the cloud-shadow reflectance spectrum is as much as 
15% below that for the conventional approach for wavelengths less than 580 nm. For 
longer wavelengths the comparison is still excellent. The area about the shadow of Cloud 
2 appears less isotropic than that about the shadow of Cloud 1, perhaps contributing to the 
larger difference. The conventional method, however, would have seriously over- 
corrected for aerosol radiance had sunglint and/or foam effects not been detected and 
separated out using the cloud-shadow method. 

Thus, although the cloud-shadow method cannot provide a correction for the entire 
image (unless the aerosol concentration and type are horizontally homogeneous), it can be 
used effectively to provide an independent check for specific locations to help identify 
problems resulting from sensor calibration or atmospheric-removal methodology. 

Appendix 

The term Ask> .L ws ( = ^yL^ - ^yL^ ) which appears in Eq. (13) cannot be derived from 
Lowtran 7. Instead, this quantity was estimated using backward Monte Carlo simulation 1 '• 
,7 ’ l8 . The modelled atmosphere consisted of 50 horizontally infinite and homogeneous 
layers overlying a flat sea surface. Optical properties of the atmosphere were taken from 
Elterman 12 , with Rayleigh and marine aerosol phase functions as in Reinersman and 
Carder". 
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The quantites of interest in these calculations were the differences in downward 
diffuse irradiances needed to calculate Photons reflected from the surface of the 

sea make some contribution to E tl (0 + ), as do photons backscattered back into the 
atmosphere from below the sea surface. But these contributions to E d (0 + ) are probably 

about the same in the shadow and neighbor regions, and would cancel each other in the 
calculation of So, for simplicity, the sea surface was modelled as a perfect 

absorber, and photons which would have impinged upon the sea surface more than once 
were neglected. 

Clouds were modelled as spheres which completely displaced the atmosphere in 
which they were imbedded. The extinction coefficient within clouds was set at 50 km 1 , 
and the single scattering albedo was set to unity 19 . Scattering within clouds was governed 
by a single-term Henyey-Greenstein phase function with the asymmetry factor set to ,85. 20 
In the results to follow. Cloud 1 and Cloud 2 represent the bottom and top cloud-shadow 
pairs of Fig. 4, respectively. Cloud 1 has a radius of 372 m and a center height of 942 m. 
Cloud 2 has a radius of 350 m and a center height of 1 194 m. These parameters were 
based on estimation of the areas of the shadows and the solar zenith angle of 50.8°. 

Photon packets were traced backward from a receiver located on the surface at a point 
of interest using a method derived from Gordon 18 . Each simulation in this work traced 10 5 
packets, each initially representing 10 10 photons, until the weight of the packet was 
diminished to less than 1 photon. Estimates of the downward diffuse irradiance 
(normalized by the extraterrestrial solar irradiance) and the average cosine of the 
downward diffuse irradiance were acquired for each sensor position modelled. 

Preliminary simulations indicated that the diffuse light field is highly variable in the 
region near and within the shadow. Figs. 8 and 9 show the normalized diffuse irradiance 
and the average cosine for diffuse irradiance at the surface in the vicinity of model Cloud I . 


13 



These plots represent results modelled at 400 nm with Elterman's standard atmosphere. 

The points in Fig. 8 are located in the vertical plane containing the sun, cloud center, and 
shadow center ( SCS plane ). The points in Fig. 9 lie on the line perpendicular to the SCS 
plane through the center of the shadow. Reference lines are included which show the 
normalized downward diffuse irradiance and average cosine for the same solar zenith angle 
and atmospheric conditions with no cloud present. 

An interesting feature of Fig. 8 is the bright spot appearing directly below and sunward 
of the cloud. This may be the result of modelling a spherical cloud, since, in this case, part 
of the brightly illuminated side of the cloud is visible from directly below. The features 
relevant to this work, however, are the depletion of diffuse irradiance at the center of the 
shadow, and the distance to which perturbations due to the presence of the cloud extend 
from the center of the shadow. Many combinations of cloud sizes, cloud heights, and 
aerosol optical thicknesses were modelled, and the same general features appeared in each 
case. 

Selection of the neighboring region for a cloud shadow involves a compromise. The 
neighboring region should be near enough to the shadow that the inherent optical 
properties of the water in both regions are the same, but the neighboring region should be 
in a location where the downward radiance distribution at the surface is not perturbed by 
the presence of the cloud. Consider the two half-planes defined by the line through the 
shadow center and perpendicular to the SCS plane. The results shown in Figs. 8 and 9 
indicate that the neighbor region should be located in the half-plane which lies further from 
the cloud, i.e. on the side of the shadow from which the illuminated face of the cloud is not 
visible. Neighbor points chosen directly on the SCS plane should be at least 5 cloud radii 
from the shadow center. Those chosen along the line perpendicular to the SCS plane may 
be as near as 3 cloud radii to the shadow center. 
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The data comprising Fig. 10 resulted from modelling Cloud 1 and Cloud 2 under the 
same conditions just described. The neighbor region was located in each case 3 cloud radii 
from the shadow center on the line through the shadow center and perpendicular to the 
SCS plane. The average cosine of the diffuse irradiance at the neighbor regions is 
indistinguishable from the clear sky ( no cloud present ) value. The average cosine at the 
shadow centers varies from the clear sky values by less than 10% in the worst case. The 
normalized downward diffuse irradiance at the neighbor points is also indistinguishable 
from the clear sky value. The normalized downward diffuse irradiance at the shadow 
centers differs between Cloud 1 and Cloud 2 due to differences in cloud sizes and heights. 

Elterman's standard atmosphere represents conditions which are much more turbid 
than those prevailing when the image of Fig. 4 was acquired. Accurate modelling of these 
very clear conditions required scaling of Elterman's aerosol profile so that t a (550 nm) = 

.0375 . The simulations for Fig. 1 0 were repeated using the scaled atmosphere, and the 
results are shown in Fig. 1 1 . Note that for the clear conditions, the difference in diffuse 
irradiance at the center of the shadows for the two clouds is negligible, and that the 
difference in diffuse irradiance between the shadow and neighbor regions is much less than 
for the turbid conditions shown in Fig. 10. 


Geometric arguments dictate that if the average cosines for downward diffuse 
irradiance were the same in the shadow and neighbor regions, then the remote sensing 
reflectances for downward diffuse irradiance in the two regions would also be equal. 
Figures 10 and 1 1 show that the average cosines are not equal, but vary by less that 1 0%. 
For the present, assume that the small differences in downward average cosines may be 
neglected. Then 


skyf^rs( shadow) sky^rs( neighbor) ^ts , 


(16) 


and 
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(17) 


AskyLws l^TS ( sky^d(0 ((neighbor) “ ((shadow) ) • 


The consequences of miscalculating Asky L ws can be understood qualitatively by examining 
its importance in Eq. (13) with respect to the term ( + Asky L ws ) . Again assuming 

thatEq. (16) holds 


AskyLws ( sol^wn AskyLws ) ( sky^dfK ((neighbor) “ ^(0 ((shadow) ) 


^ ( sol^d(d ((neighbor) sky^d(d ((neighbor) “ Ej(0 ((shadow) ( • 


(18) 


The values of the right-hand-side of Eq. ( 1 8) are shown in Fig. 12 for Cloud 1 and Cloud 2, 
each embedded in both Elterman's standard and scaled atmosphere. Note that omission of 
Asky L ws can lead to an error of up to 25% in the calculation of and hence in R re , when 

modelling clouds about the size and height used in this demonstration if they are embedded 
in a very turbid atmosphere. However, accurate estimation of Asky L ws becomes less 

important when using the cloud-shadow method in very clear conditions, where it 
contributes less than 3% in the calculation of R IS ( 700 nm ). 
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Figure Captions: 


Fig. 1 : A geometric demonstration scene of a paired cloud and shadow arrangement. 

Fig. 2: 0 (r,X), from backward Monte Carlo simulation. 

Fig. 3: A scaling relationship between L a (780 nm) and x a (780 nm) derived from 
LOWTRAN-7. 

Fig. 4: The locations of sites where cloud-shadow method was used. Shadow and 
neighbor regions are marked by boxes. 

Fig. 5: Illustration of the convergence of the cloud-shadow calibration. Solid line is "true" 
Rrs, dotted line is R' re which would be derived if sensor calibration had been in error by 

10%. Intermediate lines indicated the iterative convergence of the cloud-shadow method. 

Fig. 6: Remote-sensing reflectance spectra from the bottom site neighborhood of Figure 2 
using conventional atmospheric-correction methods and the cloud-shadow method. 

Fig. 7: Remote-sensing reflectance spectra from the top site neighborhood of Figure 2 
using conventional atmospheric-correction methods and the cloud-shadow method. 

Fig. 8: Normalized downward diffuse irradiance and average cosine for diffuse irradiance 
in the vicinity of Cloud 1 at 400 nm in Elterman's standard atmosphere. Points lie on the 
surface in the vertical plane containing the sun, cloud, and shadow center. 

Fig. 9: Normalized downward diffuse irradiance and average cosine for diffuse irradiance 
in the vicinity of Cloud 1 at 400 nm in Elterman's standard atmosphere. Points lie on the 
surface along a line through the center of the shadow and perpendicular to the vertical 
plane containing the sun, cloud, and shadow center. 
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Fig. 10: Normalized downward diffuse irradiance and average cosine at shadow center 
and neighbor regions for Cloud 1 and Cloud 2 embedded in Elterman's standard 
atmosphere. Clear sky values represent conditions with no cloud present. 

Fig. 11: Same as Figure 10, except Elterman's aerosol scaled so that t a (550 nm) = .0375 . 

Fig. 12: Values of RHS of Eq. ( 1 8) for Cloud 1 and Cloud 2 embedded in Elterman's 
standard and scaled atmospheres. 
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Figure 4. : The locations of sites where cloud-shadow method was used. Shadow and neighbor regions are 
marked by boxes. 
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Figure 5. Illustration of the convergence of the cloud-shadow calibration. Solid line is ’’true" R rs , dotted 
line is R’ rs which would be derived if sensor calibration had been in error by 10%. Intermediate lines 
indicated the iterative convergence of the cloud-shadow method. 


P. Reinersman, Applied Optics, Lasers, Photonics, and Environmental Optics. 




Cloud 1 



Figure 6. Remote-sensing reflectance spectra from the bottom site neighborhood of Figure 2 using 
conventional atmospheric-correction methods and the cloud-shadow method. 
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Figure 7. Remote-sensing reflectance spectra from the top site neighborhood of Figure 2 using 
conventional atmospheric-correction methods and the cloud-shadow method. 
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Figure 8. Normalized downward diffuse irradiance and average cosine for diffuse irradiance in the vicinity 
of Cloud 1 at 400 nm in Elterman’s standard atmosphere. Points lie on the surface in the vertical plane 
containing the sun, cloud, and shadow center. 
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Figure 9. Normalized downward diffuse irradiance and average cosine for diffuse irradiance in the vicinity 
of Cloud 1 at 400 nm in Elterman’s standard atmosphere. Points lie on the surface along a line through the 
center of the shadow and perpendicular to the vertical plane containing the sun, cloud, and shadow center. 
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Figure 10. Normalized downward diffuse irradiance and average cosine at shadow center and neighbor 
regions for Cloud 1 and Cloud 2 embedded in Elterman’s standard atmosphere. Clear sky values represent 
conditions with no cloud present. 
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Figure 11. Same as Figure 10, except Elterman’s aerosol scaled so that x a (550 nm) = .0375 . 
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Figure 12. Values of RHS of Eq. (18) for Cloud 1 and Cloud 2 embedded in Elterman’s standard and scaled 
atmospheres. 
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